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The singlet s-, d- and triplet p-wave pairing symmetries in quasi-one-dimensional organic supercon- 
ductors can be experimentally discriminated by probing the Andreev bound states at the sample 
edges. These states have the energy in the middle of the superconducting gap and manifest them- 
selves as a zero-bias peak in tunneling conductance into the corresponding edge. Their existence is 
related to the sign change of the pairing potential around the Fermi surface. We present an exact 
self-consistent solution of the edge problem showing the presence of the midgap states for p^-wave 
superconductivity. The spins of the edge state respond paramagnetically to a magnetic field parallel 
to the vector d that characterizes triplet pairing. 



PACS numbers: 74.70.Kn 73.20.-r 74.50. +r 



I. INTRODUCTION 

Quasi-one-dimensional (Q1D) conductors of the 
(TMTSF) 2 X family @ (the Bechgaard salts) are the first 
organic materials where superconductivity was discov- 
ered twenty years ago with T c w 1 K S. Abrikosov 
proposed that the superconductivity is p-wave triplet 
(3 1, because it is suppressed by nonmagnetic impurities 
[4|. Gor'kov and Jerome observed that the upper critical 
magnetic field H C 2 exceeds the Pauli paramagnetic limit, 
which is also a signature of triplet superconductivity || . 
Recent data show that H C 2 exceeds the Pauli limit by 
a factor greater than 4 0. Another signature of triplet 
pairing is that the Knight shift does not change between 
the normal and superconducting states 0. However, the 
temperature dependence of the NMR relaxation rate ||] 
and analogy with the high-temperature superconductors 
led to an alternative proposal of the d-wave symmetry || . 
It was also proposed that a singlet Q1D superconductiv- 
ity can overcome the Pauli paramagnetic limit by form- 
ing the spatially nonuniform Larkin-Ovchinnikov-Fulde- 
Ferrell state [|l0|. However, the quantitative analysis of 
the experimental data || by Lebed et al. jll],[l2| did not 
support this proposal and favored triplet pairing. The 
/-wave fl3f| was also proposed recently. So the pairing 
symmetry in the Bechgaard salts remains hotly debated. 

In this paper, we propose a phase-sensitive method to 
distinguish experimentally between the s-,p-, and <i-wave 
symmetries. We employ a relation between sign change of 
the superconducting pair potential around the Fermi sur- 
face and existence of the surface Andreev bound states, 
discovered for p-wave by Buchholtz and Zwicknagl fl^ j 
and for d-wave by Hu ]l~5| . For different superconducting 
symmetries, we determine which edges of (TMTSF^X 
must have the Andreev bound states. The energy of 
these states is in the middle of the superconducting gap, 
thus they can be observed in tunneling experiments as 
zero-bias conductance peaks |TJ,[l6|,[l7j . We also obtain 



an exact self-consistent solution of the edge problem for 
a Pa-wave Q1D superconductor by mapping it onto the 
kink soliton solution for a ID charge-density wave |fl8| . 
We show that the spins of the edge states should exhibit 
a strong paramagnetic response to a magnetic field paral- 
lel to the polarization vector d of the triplet pairing and 
propose the corresponding experiment. All calculations 
are performed at zero temperature. 



II. Q1D SUPERCONDUCTIVITY 

Classification of superconducting pairing symmetry is 
particularly simple for a ID electron gas. Its Fermi sur- 
face consists of two points ±fcp. Let us introduce the 
operators "tp" of the right (a=R) and left (a=L) mov- 
ing electrons with the momenta close to ±/cf and the 
spin cr=|,|. The Cooper pairing can be either singlet 

(CC,) a e ^' A " = l& i V v> Aa or tri P let « 

ia^ v \d ■ <r)A". Here a=L,R for ev=R,L; e CTCT < is the an- 
tisymmetric metric tensor, and <r are the Pauli matrices 
acting in the spin space; d is a unit vector of polarization 
of the triplet state. Since the fermion operators anticom- 
mute, the superconducting pair potential has either the 
same (A R = A L ) or the opposite (A R = — A L ) signs at 
the two Fermi points for the singlet or triplet pairing. 

The real (TMTSF^X materials are three-dimensional 
(3D) crystals consisting of parallel chains. In the tight- 
binding approximation, the electron energy dispersion 
(measured from the Fermi energy) can be written as |h| 



e(k) = «f(|^x| — kp) — 2tbCos(k y b) — 2t c cos(k z 



(1) 



In the right-hand side of Eq. (Q) , the first term represents 
the dispersion along the chains, linearized near the Fermi 
energy with a Fermi velocity up- The two other terms 
describe electron tunneling between the chains in the y 
and z directions with the amplitudes tb and t c . k = 
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FIG. 1. Top: (TMTSF) 2 X samples with the lines indi- 
cating ID chains. The left and right panels sketch tunnel- 
ing along the a and 6 axes. Bottom: The Fermi surface of 
(TMTSF)2X, sketched with a greatly exaggerated warping in 
the k y direction. Reflection from the edge perpendicular (par- 
allel) to the chains changes electron momentum from L to R 
(R' to R and L' to L), as shown in the left (right) panel. 



and d c is unknown. In this paper, we make the simplest 
assumption that d is a real vector pointing along the a 
axis parallel to the chains (the so-called polar state p2| ) . 
If we select the spin quantization axis along d, then the 
4x4 matrix equation (|^) decouples into two 2x2 matrix 
equations for the wave functions (u a , e aS v a ): 



—iavpd x aA a (x) 
aA a *(x) iavFd x 



n,cr 



EL 



n.a 
a,<7 



(3) 



Here the index a —1,1 takes the values ± when used as 
a coefficient. It is present in the off-diagonal term aA a 
only for triplet, but not for singlet pairing. To simplify 
equations, we omit the spin index a in Sees. Ill and IV 
and restore it in Sec. [y| Notice that Eqs. (|) and 
depend only on the ID coordinate x, because the 3D 
dispersion ([!]) in k y and k z has been absorbed into the 
definition of the 3D Fermi momenta kp. 



III. EDGE STATES 



(k x ,k v ,k z ) is the 3D electron momentum, b and c are 
the lattice spacings in the y and z directions, and U = 1. 

The Fermi surface corresponding to Eq. (0) consists of 
two disconnected sheets, sketched in Fig. ljwith a greatly 
exaggerated warping in the k y direction. In the simplest 
case, the superconducting pair potential A(k) is equal 
to a constant A Q on a given sheet a of the Fermi sur- 
face, and A R = ±A L for the singlet or triplet pairing, 
respectively. In both cases, the superconducting gap has 
no nodes on the Fermi surface. These two symmetries 
can be called s- and p x -w&ves. Other symmetries will be 
discussed at the end of the paper (see Fig. ||) . 

Electron eigenstates of the energies E n are described in 
a superconductor by the Bogolyubov-de Gennes (BdG) 
wave functions $ n = e 4r ' kp [u n _ a {x), e a >a>v° (x)), where 
«„ jCr and are the electron-like and hole-like compo- 
nents with the spins a and a', and a' =1, f is the spin 
index opposite to a' =|, J,. The 3D Fermi momenta kp 
belong to the warped Fermi surface shown in Fig. [l]. Near 
the Fermi surface a, the wave functions satisfy the lin- 
earized BdG equation pfj|,^l| 



— iavpd x (a--d)A a (x] 
(<x • d) A Q *(x) iav F d x 



E r . 



, (2) 



where avp = ±vp for a=R,L. The term (<x • d) is present 
for triplet superconductivity and absent for the singlet 
one. It operates on the spin indices of the components u 
and v. 

In general, the vector d is a function of the position on 
the Fermi surface, e.g. in 3 He-A and 3 He-B Q. The 
quantitative analysis |l2| of the experimental data 
gives the following information about the components of 
the vector d in the (TMTSF) 2 X crystal: d a ^ 0, d b = 0, 



Let us consider a system occupying the semi-infinite 
space x > with an impenetrable edge at x — 0. When 
the electron reflects from the edge specularly, its k x mo- 
mentum changes sign, whereas the other components re- 
main the same. The electron scatters from the point kp 
at the left sheet of the Fermi surface to the point kp at 
the right sheet, as shown in the left panel of Fig. [j] Thus, 
its BdG wave function \& is a superposition of the R and 
L terms: 



1 

71 



<0) 



(4) 



We have selected the minus sign in Eq. (^) so that the 
impenetrable boundary condition ty(x — 0) = gives 



M R (0)= M L (0), v R (0) = v L (0). 



(5) 



First let us use a step-function approximation for the 
pairing potential: |A Q (a;)| = Aq6(x). Then, the plane 
waves [u a (x) , v a (x)] oc e lk " x are the eigenfunctions of 
Eq. (||) with the energies E = ±y/ (vFk x ) 2 + A^. How- 
ever, the energy is real also when k x is imaginary (but 
not a combination of real and imaginary parts): k x — in 
and E = ±yA~Q — (v-pn) 2 . For k > 0, this solu- 
tion describes an electron eigenfunction localized near 
the edge at x — 0: [u a (x) , v a (x)] oc e~ KX . Because 
u a /v a = A a /(aivFH + E), the boundary condition (j|) 
can be satisfied only for p^-wave with A R = — A L , but 
not for s-wave with A R = A L . Thus, in the p x case, there 
is an edge electron state with the energy in the middle 
of the superconducting gap: E = 0, and the localization 
length is equal to the coherence length: 1/n — vp/ 'A . 

The step-function approximation does not take into 
account the BdG self-consistency condition A a (x) = 



2 



gJ2 n u Zi x ) v n* ( x ) ' wner e g is the effective coupling con- 
stant |23], and the sum is taken over all occupied states 
with E n < (at zero temperature). To solve the prob- 
lem, let us extend the wave function (Q) from the positive 
semispace x > to the full space — oo < x < oo. Let us 
define [u(x),v(x)] = [u R (x),v R (x)] and A(x) = A R (x) 
for x > 0, and [u(x),v(x)] = [u L (— x), v L (— x)] and 
A(x) = A L (— x) for x < 0. Because of the boundary 
condition (|B|), the wave function [u(x),v(x)] is contin- 
uous at x = and satisfies a single BdG equation for 
— oo < x < oo with the BdG self-consistency condition: 



-iv F d x A(x) \ fu n (x) 
A* (as) + i w F 3 x J ^y n (x) 



u„(x) 
u„(x) 



(6) 
(7) 



Eqs. (|J) and (Q) coincide with the exactly solvable 
equations describing ID charge-density wave in poly- 
acetylene 18 . The p x -wave problem, where A(x) 
changes sign: A(+oo) = — A(— oo), maps onto the kink 
soliton solution [n8|: 



A(x) = «Ao tanh(Kx) 



vq(x) J 2 cosh(Kx) \— 1/ 



E k = ±Jvlk 2 



Uk (x) 



A 2 



E k + v F k + A(x) 
Vk(x) J 2E k ^/~L~ x \E k - v F k - A(x) 



(8) 
(9) 

(10) 
(11) 



where L x is the length of the sample along the chains. 
One can check explicitly that solution (P)-(pT|) satisfies 
Eqs. (||) and (f| p3]. It also has the property of su- 
persymmetry p5| |. The localized electron state (^) with 
Eq = corresponds to the Andreev edge state in the p x - 
wave superconductor. In the s-wave case, where A(x) 
does not change sign: A(+oo) = A(— oo), the solution of 
Eqs. @ and (|) gives a uniform A(x), which does not 
have bound states. The existence of the midgap state 
in the case where A(x) changes sign is guaranteed by 
the index theorem and does not depend on the detailed 
functional form of the pair potential |2q] . 

BdG states are described by the operators i& — utp + 
v*xf)\ The expectation value of electric charge p in the 
edge states is zero: p oc |uo| 2 — |«o| 2 = 0. For a ID p x - 
wave superconductor with only one species of spin, Eqs. 
(Q) and @ imply that the two edge states at the oppo- 
site ends are described by the Majorana operators of the 
opposite parity: iff' = i^P |27|. There was a proposal 
to use such edge Majorana fermions for quantum com- 
puting p7| . However, in a Q1D p^-wave superconductor, 
the midgap states with different momenta k|| parallel to 
the edge and spins a form a degenerate continuum with 



IV. TUNNELING 

Let us consider electron tunneling between the super- 
conducting (TMTSF)2X and a normal metallic tip. The 
tunneling junction can be modeled as two semi-infinite 
regions, normal (N) and superconducting (S), with a flat 
interface between them. Following Refs. |I7J,^9|,[30|] , we 
solve the BdG equations in the ballistic regime assuming 
specular reflection and the translational invariance paral- 
lel to the interface J3lj . To make the problem analytically 
tractable, we use the step-function approximation for the 
pair potential. At the interface, we impose the boundary 
conditions * N = * s and u n *n = #s*S + 2i7Y* N @, 
where £>n,s are the components of the velocity operators 
perpendicular to the interface in metal and superconduc- 
tor, and Ti. is the strength of the interface barrier. From 
the solution of the BdG equations, we find the probabil- 
ities S(£',k||) and A{E, kn) of the normal and Andreev 
|Uf reflections as functions of the electron energy E and 
momentum k|| parallel to the interface. They determine 
the dimensionless conductance G = l + A — Bin the 
formula for the electric current through the contact [E9J : 



I 



2eS 
~h~ 



d 2 k v dE 



(2tt) 2 



[f(E-eV)- f{E)]G{EM). (12) 



Here S is the contact area, V the bias voltage, f(E) the 
Fermi function, e the electron charge, and h the Planck 
constant. It follows from Eq. ( |l2"| ) that the differential 
conductance of the contact at zero temperature, 



G{V) 



dl 
dV 



2e 2 S f d 2 fc|. 



{2-Kf 



G{eVM 



(13) 



is proportional to the average over k|| of the dimension- 
less conductance G |53| . The latter is determined by the 
transmission coefficient T at a given k|| |l6| , |32| ]: 



G+=T 



1 + T\T\ 2 + {T -1)\T\ 4 
|l±(T-l)r 2 | 2 



(14) 



where 



r , F ,_ / [E-sg n(E)y/E* -Al]/A , \E\>A , 
\{E-i^Aj^, W)/A , \E\<A , 



T = 4v N v s /[(v N + v s ) 2 + m 2 }. 



(16) 



The ± sign in Eq. (|TJ) is the relative sign of the pair 
potentials for the two branches of BdG quasiparticles in- 
volved in tunneling. For tunneling along the chains, the 
two branches correspond to the points L and R in the left 
panel of Fig. [l], and the sign in Eq. ( |l4] ) is sgn(A R A L ): 
+ for s-wave and — for p^-wave. Averaging in Eq. ( |T3| ) 
is performed taking into account that vjy & n d v§ in Eq. 
( p6| ) and Ao in Eq. ( |l5| ) may depend on k|| . 

As follows from Eq. (|l4|), G+ and G_ coincide for a 
fully transparent interface (T = 1): G+ = G- = 1 + 
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FIG. 2. Dimensionless conductances G_ (left panel) and 
G+ (right panel) given by Eq. ( |l4| ) are plotted versus energy 
E for the transmission coefficients T — 0.5 and 0.25. G+ 
and G- correspond to the cases where the superconducting 
pairing potential has the same or the opposite signs for the 
two branches of BdG quasiparticles involved in tunneling (the 
points L and R, R' and R, L' and L in Fig. |l|). 



|T| 2 |||. However, typically T < 1, both because of 
the barrier potential H. and the mismatch of the normal 
Fermi velocities vn ^ v$ in metal and superconductor in 
Eq. (|l|) HI). At low interface transparency T < 1, G_ 
and G+ behave as shown in the left and right panels of 
Fig. H for T = 0.5 and T — 0.25. Inside the energy gap, 
where |T| = 1, G-(E) has a Lorentzian shape with the 
maximum of 2 at E = 0, the width proportional to T, 
and the minimum proportional to T 2 at \E\ = Aq: 



G-(E) 



T 2 /2{\-T) 



(£/A ) 2 + T 2 /4(l-T)' 



|£|<A . (17) 



G + (E) shows the opposite behavior: a minimum propor- 
tional to T 2 at E = and the maxima of 2 at \E\ = Aq: 



G+(E) 



T 2 /2{l-T) 



l-(^/A ) 2 +T 2 /4(l-T)- 



\E\ < A . 



Both G+ and G_ approach the normal-state conductance 
T at \E\ > A . 

The zero-bias conductance peak (ZBCP), shown in the 
left panel of Fig. ||, is a manifestation of the midgap An- 
dreev bound states. They exist at those edges where 
momentum reflection from the edge connects the points 
on the Fermi surface with opposite signs of the super- 
conducting pair potential. As shown in the left and right 
panels of Fig. |l|, reflection from the edge perpendicular to 
the chains connects L to R, and reflection from the edge 
parallel to the chains connects R' to R and L' to L. By 
comparing the signs of the pair potential at these points 
for the superconducting symmetries |l2| listed in Table 
| and sketched in Fig. ||, we determine whether ZBCP 
must be present in tunneling into those edges. Compar- 
ison of Table || with the experiment should uncover the 
superconducting symmetry of (TMTSF) 2 X. 



TABLE I. Presence (Yes) or absence (No) of a zero-bias 
conductance peak in electron tunneling along the a and b axes 
(see the top left and right panels in Fig. |l|) for different sym- 
metries of the superconducting pairing potential A(k) p^ , p5[ 



Symmetry 


A(k) 


a-axis ZBCP 


6-axis ZBCP 


s 


const 


No 


No 


Px 


sin(fc^a) 


Yes 


No 


Pv 


sin(fcyfe) 


No 


Yes 


^x 2 —y 2 


cos(k y b) 


No 


No 




sva(k x a) sm(k y b) 


Yes 


Yes 



V. SPIN RESPONSE 

Now let us discuss the spin response of the edge states 
in a triplet superconductor subject to an external mag- 
netic field H. (Here we do not consider orbital effects 
of the magnetic field, such as the Meissner effect.) In 
this case, the matrix in Eq. (||) should be replaced by the 
following matrix: 



-iavpd x — fiB (H • tr) (<x • d) A a (x) 
(<r ■ d) A a * (x) iav F d x - fj, B (H • tr) 



(18) 



where /is is the Bohr magneton, and the electron (/-factor 
is 2. 

If H||d, then, by selecting the spin quantization axis 
along d, the 4x4 matrix equation (|lj) is decoupled into 
two 2x2 matrix equations similar to Eq. (||) for the wave 
functions (u a , e aS v a ) with the matrix 



-iavpdx — u^bH aA a (x) 
aA a *(x) iavpdx — <j[IbH 



(19) 



For a given spin projection <r, the magnetic field H enters 
Eq. ( |l9| ) as a unity matrix and simply shifts the spectrum 
(0) and ( |Io| ) by —u^lbH. Thus, the energies of the up 
and down spin states become split by =F/isi?, includ- 
ing the midgap state: Eq = ^f/isH. Because this state is 
half-filled, the edge states with the spin parallel (antipar- 
allel) to the magnetic field become completely occupied 
(empty). This generates spin h/2 and magnetic moment 
/i^ at the end of each chain. Such a giant magnetic mo- 
ment was predicted by Hu |l5],[36| for the edge states in 
a singlet d-wave superconductor. In a triplet supercon- 
ductor, the effect is similar, but anisotropic. 
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FIG. 3. Different symmetries of the pairing potential A(k) 
in a Q1D superconductor. The solid and dotted lines repre- 
sent the portions of the Fermi surface with the opposite signs 
of the pairing potential. 
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Indeed, suppose now that Hid. In this case, it is 
convenient to select the spin-quantization axis z along H 
and the x axis along d. Then Eq. @ separates into 
two 2x2 matrix equations similar to Eq. (|J) for the wave 
functions (u a , e Sa v' 7 ) with the matrix 



d 



(20) 



-iav-pdx — u/ibH A a (x) 

A a *(x) iavpd x + u^bH 

The magnetic field H can be eliminated from Eq. (|2C 
by adjusting the Fermi momenta for the up and down 
spin states: kp,a = kp + <j^bH /vp- Thus, the energy 
spectrum of the system remains the same as in Eqs. (^) 
and (|l0|). Particularly, the energy of the midgap state 
does not split: Eq = 0, thus no unbalanced spin and 
magnetic moment are generated on the edge. 

We see that the edge spin response of a triplet su- 
perconductor is opposite to its bulk spin response. It 
is well known p2| that the bulk spin susceptibility for 
H _L d is the same as in the normal state, whereas for 
H||d it vanishes at zero temperature. For the edge states, 
the spin response vanishes for Hid and is paramag- 
netic for H||d. Nominally, the edge spin susceptibility 
is infinite, because, formally, an infinitesimal magnetic 
field can completely polarize the edge spins. We can 
only estimate the maximal magnetic moment, which is 
(i B = 9.3 x 1(T 24 A m 2 per chain or fi B /bc = 9x 1CT 6 
fiA per unit area of the edge, where b — 0.77 nm and 
c = 1.35 nm @. 

The generation of paramagnetic moment by the edge 
states for H||d could be observed experimentally by mea- 
suring magnetic susceptibility with a coil as shown in Fig. 
[|, where we assume that d is directed along the chains. 
In the bulk, far from the edges, the susceptibility should 
be diamagnetic, because of the orbital Meissner effect j37j 
and vanishing spin bulk susceptibility for H||d. However, 
when the coil is moved toward the sample end, the sus- 
ceptibility should change sign and become paramagnetic 
because of the edge states. They are localized within the 
coherence length £ = Tivf / Aq = 0.6 /xm, where we used 
A = 0.22 mcV, v F = t a a/V2h = 190 km/s, t a = 0.25 
eV, and a = 0.73 nm fisfl . The effect should depend on 
the coil orientation relative to the vector d. Therefore, 
this experiment could confirm the existence of the edge 
states and the pairing symmetry in the (TMTSF^X su- 
perconductors. 




H 



FIG. 4. Schematic experimental setup to measure magnetic 
susceptibility of the edge states localized at the ends of the 
chains. 



themselves as zero-bias peaks in tunneling conductance. 
Thus, the pairing symmetry of the Q1D superconductors 
can be determined by tunneling into the edges perpen- 
dicular and parallel to the chains. The spins of the edge 
state respond paramagnetically to a magnetic field par- 
allel to the vector d that characterizes triplet pairing, 
generating the magnetic moment fis per chain. 

The (TMTSF)2X materials are expected to have elec- 
tron edge states also in the magnetic-field-induced spin- 
density- wave phase (FISDW), which exhibits the quan- 
tum Hall effect p8| . Those states are chiral and have 
dispersion inside the energy gap. The midgap states dis- 
cussed in the present paper also acquire chiral dispersion 
due to the orbital effect of a magnetic field, similar to 
cuprates |39|. Another interesting example of the edge 
states in a ID electron gas was studied theoretically in 
Ref. ]4C| ] for the Luther-Emery fermions in the bosonized 
representation. 

K.S., H.J.K., and V.M.Y. were supported by the 
Packard Foundation and NSF Grant No. DMR-9815094; 
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Note added in proofs. A more subtle consideration of 
the results of Sec. V shows that the spin of an edge state 
is actually fractional and is equal to h/A (per each end 
of each chain), not h/2. Correspondingly, the magnetic 
moment is /is/2, and the numerical estimates given in 
Sec. V should be multiplied by an additional factor 1/2. 
Derivation of these results will be given in a separate 
paper. The authors are grateful to A. Yu. Kitaev and D. 
A. Ivanov for very illuminating discussions of this subject 
(see D. A. Ivanov, |cond-mat/9911147| ). 



VI. CONCLUSIONS 

We have constructed an exact analytical self-consistent 
solution of the edge problem for ap^-wave Q1D supercon- 
ductor by mapping it onto the kink soliton solution for a 
ID charge-density wave. The edge electron midgap states 
exist when the pairing potential has opposite signs at 
the different parts of the Fermi surface connected by mo- 
mentum reflection from the edge. These states manifest 
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